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We study the Langevin dynamics of a ferromagnetic Ginzburg-Landau Hamiltonian with a competing long- 
range repulsive term in the presence of an external magnetic field. The model is analytically solved within 
the self consistent Hartree approximation for two different initial conditions: disordered or zero field cooled 
(ZFC), and fully magnetized or field cooled (FC). To test the predictions of the approximation we develop a 
suitable numerical scheme to ensure the isotropic nature of the interactions. Both the analytical approach and 
the numerical simulations of two-dimensional finite systems confirm a simple aging scenario at zero temperature 
and zero field. At zero temperature a critical field h c is found below which the initial conditions are relevant for 
the long time dynamics of the system. For h < h c a logarithmic growth of modulated domains is found in the 
numerical simulations but this behavior is not captured by the analytical approach which predicts a t 1//2 growth 
law at T = 0. 



I. INTRODUCTION 

It is now well understood that modulated structures can 
arise from underlying competing interactions acting upon dif- 
ferent scales QD. In magnetic systems, e.g. in ultrathin mag- 
netic films with perpendicular anisotropy HHH], these phases 
are present due to an exchange (ferromagnetic) short range 
interaction competing with a weaker but long range (antifer- 
romagnetic) dipolar interaction. This competition gives rise to 
striped and bubbled patterns that break the translational and/or 
rotational symmetry of the system. 

This phenomenology is shared with a large number of sys- 
tems. Type-II superconductors |4], doped Mott insulators 
Bh, quantum Hall systems |0] and diblock copolymers 10, 0] 
are other examples of systems with competing interactions 
which have symmetry breaking transitions to inhomogeneous 
phases. Since domains in these systems often are of meso- 
scopic scales, it is possible to perform a general analysis 
which neglects the microscopic details of each specific system 
while capturing its general properties. A pioneering attempt 
in this direction was done by Brazovskii EJ], who showed that 
systems in which the spectrum of fluctuations attain a mini- 
mum in a shell at a non-zero wave vector, undergo a first order 
phase transition to a modulated state driven by fluctuations, in 
contrast to the second order transition predicted by the mean 
field theory. Moreover, the presence of competing interac- 
tions frequently lead to a proliferation of metastable states at 
low temperatures that may dominate the dynamical behavior 
of the system. For the three dimensional case, recent exper- 
imental results in colloidal systems [10] report the presence 
of glassy states, in agreement with theoretical predictions for 
a model with screened electrostatic repulsion 111 ill and for a 
Coulomb frustrated ferromagnet Il^l - [l4ll . 

However, it is still far from clear which are the relevant 
aspects that control the long time dynamics of these systems 



and, more generally, what kind of dynamics should one ex- 
pect. In an early work for a model with ferromagnetic and 
dipolar interactions, Roland and Desai H 1 311 studied the dy- 
namics in the short-time regime, showing that the order pa- 
rameter evolves to modulated structures. Almost contempo- 
raneously, Elder and Grant [11611 obtained a scaling relation 
for the dynamics of phase ordering of the modulated phases, 
characterized by a growth law L ex t x l z with a dynamical 
exponent z = 2. After that, many works were devoted to ver- 
ify such scaling relation, solving numerically the dynamics of 
the Swift-Hohenberg model |17L Indeed most of these 
numerical simulations fl9l423ll report a late-time regime that 
satisfies scaling, but with a growth exponent near z = 5 for 
zero temperature quenches and z = 4 for low finite temper- 
atures, obtained mainly monitoring the structure factor and 
related quantities. In 112411 an exponent z = 3 was observed 
in zero temperature quenches which was associated to grain 
boundary motion. For quenches at zero temperature, a dy- 
namic crossover to a frozen state wa s rep orted 1 25], and asso- 
ciated with grain boundary pinning 11261 . On the other hand, 
recent experiments in diblock copolymers Il27ll (that are as- 
sumed to be good model systems for studying phase order- 
ing of modulated structures), show that the late time dynam- 
ics of the stripe ordering satisfies scaling with z = 4, driven 
mainly by the annihilation involving tripoles and quadrupoles 
of disclinations. For a bubble forming system Harrison et al. 
1 28] found the same growth law, but this time dominated by 
the collapse of smaller grains which reside on the boundary 
of two larger grains. More recently, Gomez et al. Ii2911 found 
a logarithmic growth L(t) ~ hit of domains in simulations 
of an hexagonal diblock copolymer model. The logarithmic 
growth is attributed to the Lifshitz mechanism of pinning of 
triple points at grain boundaries. Here we report logarithmic 
growth of different characteristic length scales in simulations 
of a stripe forming system. Also, Gleiser et al. II 3 Oil found 
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a crossover between an apparently logarithmic relaxation to 
at 1 / 2 law in the coarsening dynamics of a dipolar frustrated 
Ising model. So, despite a lot of numerical and experimen- 
tal effort the question of dynamical universality of modulated 
systems remains elusive. 

Recently, Mulet and Stariolo Hi] solved the long time 
Langevin dynamics of a model with competing interactions in 
the Hartree (self-consistent field) approximation. For the case 
of zero external field they showed that the dynamics can be 
separated in two times scales. In the first one modulated struc- 
tures emerge, and the behavior resembles the formation of typ- 
ical ferromagnetic domains [ 32] . Once these modulated struc- 
tures are formed the dynamics changes qualitatively, becom- 
ing independent of the system dimension and the temperature 
is a relevant variable. For T > the system exhibits inter- 
rupted aging and a standard paramagnetic (exponential) relax- 
ation for large times. At low temperatures domains of stripes 
of finite size, characterized by a finite correlation length, are 
formed. At T = the correlation length diverges and stripes 
order sets in with a growth law L ex t 1 ! 2 . Our simulations 
show that these results are correct in the time scale while the 
very slow dynamics of topological defects does not play a sig- 
nificant role. 

In this work we extend the results of ll3lll by solving, within 
the Hartree approximation, the same model in the presence of 
an external magnetic field and for two different initial condi- 
tions. The initial conditions were chosen to mimic standard 
experimental protocols: in the Zero Field Cooling (ZFC) ex- 
periment, the system is initialized in a highly disordered state 
(high temperature and zero external field). After some waiting 
time the temperature is suddenly lowered towards the mod- 
ulated phase and the field is turned on. In the Field Cool- 
ing (FC) protocol, the sample is first fully magnetized in the 
presence of an external field, and then the temperature is de- 
creased keeping the field on. To check the validity of our 
approximations, we develop a reliable numerical scheme to 
reach the asymptotic dynamical behavior, while ensuring the 
isotropic nature of the model interactions. We compare the 
analytical results with extensive numerical simulations of a 
two-dimensional magnetic system with dipolar interactions. 

The rest of the paper is organized as follows. In section HI 
we present the model and the Hartree approximation to solve 
its dynamics. We present analytical results for the autocorre- 
lations and spatial correlation functions for zero temperature 
(T = 0) and finite external field (H). Section [Till is devoted 
to the numerical simulations. It begins with an analysis of 
the fluctuation spectrum of the model, which is essential to 
compare the simulations of finite systems and the analytical 
results. Then we present the numerical scheme used and com- 
pare the simulations with the analytical predictions. Finally in 
section [IV] we summarize and present the conclusions of our 
work. The details of the calculations are left for the appen- 
dices. 
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(j) A (x) — h(x)(j)(x) 



^ J d d x d d x' 0(x) J(x , x') 0(f) (1) 



where r < 0, u > and J(x,x') = J(\x — x'\) represents 
a repulsive, isotropic, competing interaction. h(x) is an ex- 
ternal field and the parameter 5 measures the relative intensity 
between the attractive and repulsive parts of the Hamiltonian. 
In the limit 5 —> oc one recovers the ferromagnetic O(N) 
model (for N = 1) 132151. 

Following the general formalism presented in Ref. Bill , the 
Langevin dynamics for Hamiltonian (Q]) takes the form: 



1 d(j){x,t) 

r m 



V 2 0(f , t) - r</)(x, t) -u <j) 3 (x, t) + h(x, t) 
- 5 j d d x f J(x,x')0(x f ,t) + (2) 



where rj(x, t) represents a Gaussian noise that models the cou- 
pling to a heat bath. 

Within the self-consistent (Hartree) approximation the non- 
linear term 3 is substituted by 3(0 2 (x, i))(/)(x, t) where the 
average is performed over the initial conditions and noise re- 
alizations. It is worth to note that, within the approximation, 
(0 2 (x, t)) is spatially homogeneous, which implies a restric- 
tion in the possible ensembles of initial conditions taken to 
perform the averages. 

We have studied two kinds of initial conditions: disordered 
(Zero Field Cooled) and ferromagnetic (Field Cooled). Both 
can be written in the general form 



(0(f,O)) = 0o 
(0(£,O)0(£',O)) = v{x-x') 



(3) 



In this way, the disordered condition (h = and T » 0) 
is obtained with 0o = and v(x — x') = AS(x — x'), while 
the ferromagnetic initial condition (h » 0) is represented by 
0o > and v(x — x') = 0q. 

After applying Fourier transforms, the equation of motion 
d21) in the Hartree approximation reads: 



at 



-[A(k) + I(t)} 0(/c, t) + 7)(fe, t) + ft(fc, t), (4) 



with 



I(t) =r o +#<0 2 (£,£)), 



A(k) = k 2 + - J(k) - a , 



(5) 



(6) 



II. ANALYTICAL APPROACH 



We consider an effective Ginzburg-Landau Hamiltonian 



where g = 3u, do = + 1/5 J(fco) and ro = r + clq. 
ko is the wave vector minimizing the spectrum of fluctuations 
A(k). In the case of isotropic competing interactions it is well 
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established that this quantity has a minimum on a spherical 
shell of radius fco in reciprocal space |9]. Thus, definitions © 
and © ensure that A(ko) = 0. 

The general solution of the dynamical equations dU) can be 
written as: 



4>(k,t) = 0(fc, 0) R(k, t, 0) + / R(k,t,t f )fj(k,t f )dt f 



R(k,t,t')h(k, t')dt' 



where R{k, t,t') = Y^y e A ^ k ^ f *') is the response function 
and Y(t) — e^o dt'i(t') # -phe Gaussian thermal noise rj(x,t) 
satisfies 

<*KM)> = o, 
(fj(k,t)fj(k f ,t f )) = 2rT(27r) d 6(k + k')6(t-t') (7) 

and a stationary homogeneous external field (h(x, t) = h) is 
defined by: 

(ft(ife)) = /i(27r) d (5(fc), 
(h(k)h(k')) = h 2 (2it) 2d 5(k)5(k'). (8) 

A complete solution amounts to the determination of the 
response function, or equivalently, Y(t). Defining K(t) = 
Y 2 (t) and following standard procedures 13 ill , it is possible to 
reduce the problem to the solution of the following differential 
equation: 
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dt 
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V(t) contains the information about the initial conditions, 
and the last two terms in © reflect the presence of the external 
field. This external field is associated to strongly non-linear 
contributions in K(t) that prevents us from using a standard 
Laplace transformation [31] to solve equation (|9]). To deal 



with these non-linearities we propose an appropriate ansatz 
for the long time behavior of K(t) and then, keeping consis- 
tency, follow usual techniques for solving differential equa- 
tions. We introduce the ansatz: 



f K^^e-^-^dt' = (K^ 2 (t) 
Jo 



(12) 



and solve for K(t), determining £ self-consistently. With this 
ansatz, the problem has been reduced to find a solution for a 
linear differential equation in K(t), a function that encloses 
conditions ©, © and ([5]). Once K(t) is known, observables 
like response functions, two times autocorrelations, and spa- 
tial correlation functions can be readily obtained. Spatial cor- 
relations are given by: 



C(x,x'^t) = C(R,t) = J ^ 



d d k' 
(27r) d 



(2n) d 

iik-x+k'-x') (13) 



where R = \x — x'\, while two-times autocorrelations are 
defined as: 



C(f,£,M') = C{t,t') = J 
d d k' 



d d k 
(2ir) d 
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^ k + k ')^ (0(k,t)4)(k f ,t f )). (14) 



Below we present and discuss the analytical predictions for 
some of these observables and leave to AppendixlAlthe details 
of the calculations. 



A. Observables at T = 

At zero temperature and in both the disordered (zero field 
cooled) and ferromagnetic (field cooled) quenches, we have 
found a critical field which divides the analysis in three re- 
gions h < h c , h = h c and h > h c , with qualitatively different 
long time dynamics ( the critical filed is defined in Appendix 
fAl). Starting from the ferromagnetic configuration we show 
that, despite the existence of this critical field, the decay of 
the correlation functions are essentially exponential, although 
with different exponential or constant rates. However, if the 
quench is done starting from a disordered high temperature 
phase, we find important functionality changes in the two- 
time and the spatial correlation functions as a function of h. 



1. Zero Field Cooling 

In table U we show the predicted long time behavior of the 
spatial and two-times autocorrelation functions in ZFC exper- 
iments. 

Coo is a constant which is different in each expression and 
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h C(M') C(R,t) 

h<h c Coo + f(R) + Coo 

(t + t')2 

h = h c Coo + r i r + Coo 

tegptpgi 

h>h c Coo+ e " Ht+ y f(R)e~ xt /Vt + Co G 
I I 

TABLE I: Long time solutions for C(£, *') and C(R, t) at T = for 
different values of the applied field in ZFC case. 



where J n (x) is a Bessel function of the first kind. 

As can be seen, for h < h c autocorrelations show aging 
behavior. In particular, the asymptotic (t^> t') two-times cor- 
relation function has the form C(t,t') oc (L(t) / L(t')) u that 
characterizes a coarsening dynamics caused by the competi- 
tion of domains of modulated structures HH). The typical size 
of these domains grows as L(t) ex t 1 / 2 , as can be confirmed 
from the exponential term of f(R) in the h < h c spatial cor- 
relations. On the other hand, for h > h c the autocorrelations 
relax exponentially fast to the asymptotic value . Also, it is 
worth to note the independence of the two-times autocorrela- 
tions with the dimensionality of the system. 

The spatial correlations at zero temperature show that only 
for h < h c the system establishes modulated structures at in- 
finite time. 



A(/i c ) = defines the critical field. Above this value of h the 
decay rate depends only on A(0) and measures, essentially, 
how fast the system relaxes to the asymptotic ferromagnetic 
state. 



III. NUMERICAL SIMULATIONS 

With the aim to test the predictions of the Hartree approxi- 
mation we solved numerically the dynamical equations Q for 
the particular case of a two dimensional model of magnetic 
film with competing exchange and dipolar interactions in the 
limit of strong perpendicular anisotropy, in which case the in- 
teractions are perfectly isotropic 113 511 . The continuum nature 
of the model and the isotropy of interactions pose some seri- 
ous challenges to the accuracy and interpretation of numerical 
solutions. We devote this section to discuss our strategy in or- 
der to obtain meaningful results. We start the section with a 
discussion of the spectrum of fluctuations of the model, how it 
is affected by finite size and discretization of space and time, 
that are inevitable when solving differential equations numer- 
ically, and propose some strategies to minimize these effects. 
After that, we describe the actual numerical scheme used to 
solve the equations. Finally, the numerical results are dis- 
cussed and compared with the analytic ones from previous 
sections. 



2. Field Cooling 

Because the initial condition is homogeneous ((j)(x) is con- 
stant), the correlations in this case are independent of R. 

As already mentioned, in this case K(t) is essentially expo- 
nential for any value of h ( see Appendix IA ll for more details). 
Moreover, the dominant exponential terms depend strongly on 
h. In particular, for h = one finds that the form of C(t, t') 
depends on the ratio between ro and -4(0), as it is shown in 
tableHH 



h 


ro 


C(R,t,t') 


h = 


r > -A(0) 


( -7 oo6 -(^o+A(o))(t+t / ) 


h = 


r < -A(0) 


Coo + e 2(ro+A(0))t 


h > 




+ e -i(A+A(0))(t+t / ) 



TABLE II: Long time solutions for C(R, t, t') and C(R, t) at zero 
temperature for FC protocols. 

The fact that auto-correlations decay to zero or non-zero 
constants according to the values of ro and A(0) can be un- 
derstood in terms of the stability of the component <j)(k = 0). 
For ro > —A(0) the k = mode is stable and the system 
evolves to a non-magnetized state. When ro < — A(0) the 
instability of the zero mode drives the system to a non-zero 
magnetization. 

For h > 0, the decay rate depends on a parameter A that 
encloses the dependency on h (see Eq. (|A19l )). The point 



A. The spectrum of fluctuations 

We discretize the system in a square mesh of L x L sites. 
In the discretized reciprocal space, the function f(t) defined 
in (fTTT) for the continuum system takes the form: 

k 

where 

fc= U^j *'J = -2 +1 '-'2' (17) 

a is the linear size of the mesh and A (k) is given by ©. From 
(IT6b we should note that only for L ^> 1 and a sufficiently 
small value of the mesh a the discrete version of f(t) can be 
approximated by the continuum one. 

At this point it is worth to note that, even in the continuum 
formulation, the functionality f(t) ~ 1/t 1 ^ 2 is only valid for 
t ^> r c , in which the function under the integral (fTTT) is sharp 
enough (see Appendix lAl and ll3lll ). Moreover, for sufficiently 
long times, the discrete nature of A(k) provokes deviations of 
f{t) from the expected power law behavior. The characteristic 
time in which this deviation happens should be of the order of 
the inverse of the difference between A(ko) and A(k), with k 
a nearest neighbor vector of ko . 

Up to now we have not specified the nature of the compet- 
ing interaction J(x, x') in (Q}. In ultrathin ferromagnetic films 
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with strong perpendicular anisotropy, the competing interac- 
tions are the short range exchange (represented in the contin- 
uum model by the square gradient term in (Q])) and the long 
ranged dipolar interaction. In the limit of strong perpendicu- 
lar anisotropy, the dipolar moments point (except at domain 
walls) out of the plane of the sample, and the magnetization 
can be modeled by a continuum scalar field <j)(x). In this limit 
the dipolar term takes the form J(x,x') oc l/\x — x f \ 3 . Note 
that in this limit the interactions are perfectly isotropic and 
so will be the spectrum of Gaussian fluctuations. The long 
range nature of dipolar interactions is commonly taken into 
account by means of the so-called Ewald Summation Tech- 
nique [ 36] . Unfortunately, the use of the J(k) given by the 
Ewald sums provokes a growing anisotropy in the fluctuation 
spectrum as the wave vector k approaches the value n/a. This 
makes the minima of the fluctuation spectrum to remain in a 
finite number of directions. In this way, after certain time, the 
parabolic nature of the isolated minima provokes a crossover 
in the exponent of the power law in f(t) from —1/2 to a di- 
mensional dependent value —dj 2, as it is shown in figure [T] 
for d — 2. This new exponent coincides with the one obtained 
for the purely ferromagnetic system [HH [37|], where f(t) is 
dominated by the single minimum of the fluctuation spectrum 
at k = 0. 
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FIG. 2: Behavior of f(t) for a system of L = 1960 using the J(k) 
given by expression dlSl . The behavior of K{t) for h = is also 
shown. The deviation from power law for large times is a finite size 
effect. 



Figure [2] shows the behaviors of f(t) and K(t) using the 
J(k) given by dT8b . The range of times in which f(t) follows 
a power law with the predicted exponent is determined only by 
r c and an upper bound dependent on the size of the system. 
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FIG. 1: Behavior of /(£) for dipolar frustrated ferromagnetic model 
of size L — 1200, using the J(k) given by the Ewald sums. For 
sufficiently large sizes a crossover is observed in the power law ex- 
ponent from — 1 / 2 to — 1 . 

In general, this long time behavior must appear in any sys- 
tem characterized by several isolated minima in the fluctua- 
tion spectrum. In this sense, it is important to be aware of this 
subtle anisotropic effect introduced by the Ewald technique. 

In order to avoid these undesired effects, a possibility is to 
replace J (k) by its expansion for small k up to second order, 
as has been already made in literature [38]. However, a more 
justified approach in the context of magnetic thin films is to 
consider the two-dimensional Fourier transform of the kernel 
of the actual dipolar interaction 1/r 3 . That is 



2tt 

J(k) = —2irk H F 

a 



1 

2' V2' 



1 



:k 2 a 2 



(18) 



where F is the confluent hypergeometric function and a is a 
short distance cutoff. 



B. Numerical scheme 

In order to build a numerical scheme it is useful to express 
© in a dimensionless form. A natural choice for the char- 
acteristic length is the wave length corresponding to the min- 
imum of the spectrum of fluctuations Ao = 27r/fco- Taking 
J(x, x') as the dipolar interaction with a short distance cut-off 
a described in the last subsection, this characteristic length is 
given by Ao = 25 (1 + ira/25). Defining also a characteristic 
time to — Aq/T and a characteristic field value (j) s = y/r/u, 
we can define new reduced variables as 



x _> t 5 a (j) 
> x, > t, yd, > a, — -> 

Ao to Aq Aq <p s 



(19) 



Equation can be rewritten as: 

a ^(f'*) = V 2 0(£,fO-i / ' d 2 x'J(\x-x'\,a)4>(x',t) 

at 5 j 

+ p[(p(x,t)-(p 3 (x,t)]+h+V2Tr](x,t) (20) 
where we have defined the dimensionless variables 



(21) 



and the thermal noise has been normalized in such a way that 

(r](x, t)r](x',t')) = 5(t — t')5(x — x'). So, we can write in 
the reciprocal space 



dct>(k, t) 
dt 



k 2 + -J(k) - p 



h - (3(p 6 (x,t) + V2Tt](x,t) 



i F 



(22) 
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where ] ? means the k component of the corresponding Fourier 
transform. 

We discretize the space with a mesh size a = 1/M. Since 
Ao is our characteristic length, this implies that the wave- 
length of the modulated structures consists in approximately 
M sites of the simulation system. Now, if we want the system 
to be able to show N modulated structures, we must simulate 
a square lattice with a linear size of L = MN sites. Thus, 
we construct a scheme in which M and N are fixed at a large 
enough value to ensure realistic simulations for the contin- 
uum model. In this sense, both quantities play a well-defined 
role. Increasing M is associated with an improvement of the 
domain wall dynamics as the magnetic structures can extend 
over several mesh sites. This results in a better definition of 
the time r c in which f(t) reaches its long times behavior. On 
the other hand, increasing N allows us to see a larger fraction 
of the (ideally infinite) system and is related to the growth of 
the maximum time for which f(t) can still be approximated 
to its long time functionality. 

In order to further simplify the spatially discretized dimen- 
sionless equation we now repeat the transformations in equa- 
tion (|T9l) for both space and time variables using a in place of 
the characteristic length (ka — >> k and t/a 2 — » t), and redefine 
the strength of the local potential as a 2 /3 — >> f3. Finally we 
consider a standard first-order semi-implicit spectral integra- 
tion scheme, with a time step dt, to construct the following 
recurrence relation: 



alytical calculations as 



r = A(k ) 
A(0) = A(0)-A(k ) 
9 = 3/3 



(23) 



where A(k) is the numerical fluctuation spectrum. In this way 
the critical field can be written as 



hc = ^^^[A(0)-A(k )] (24) 

that, under our particular set of parameters takes the value 
h c = 4.53. 



C. Simulation results 

In our simulations, the system is initialized in a disordered 
(homogeneous) configuration corresponding to the ZFC (FC) 
experiment. Suddenly, a subcritical temperature is fixed and 
the system is let to evolve a time t w , simulating a quench. 
Then, autocorrelations are computed as function of the time 
At elapsed after t w . 



1. Field Cooling 



4>(k,t + dt) = 



l + k 2 dt 



dt 



l + dt((3-aJ(k)/5) 



1 + k 2 dt 

where the noise term 77 is a random Gaussian number with unit 
variance. Here we are taking advantage of the isotropic form 
of the laplacian in the inverse space, where the adimensional 
wave vector has the form k = (2tt i/L, 2n j/L) . 



1. Choosing simulation parameters 

Simulations were carried out for a system of M = 35 and 
N = 56, that is L = 1960, with the parameters a = 0.04, 
/? = 0.265. From the chosen form of the dipolar interaction 
we have S fixed as S = (1 — ira)/2. The value of M was cho- 
sen to be large enough as to guaranty the convergence of r c to 
the continuum system value. To have a good estimate of the 
simulation time needed to reach the asymptotic regime of the 
model, we first solved numerically equation © (see figure©. 
One can see that the time at which the asymptotic function- 
ality of K(t) sets in has a lower bound different from that of 
f(t), this bound is dependent on the whole set of parameters 
and is greater than r c . 

Comparing our numerical scheme with equation © it is 
possible to express the parameters ro, A(0) and g used in an- 



We present first, the dynamical behavior of the system after 
a quench from a fully magnetized state for T = and h < h c . 
In particular, we focus our attention of the differences between 
the cases ro > —^4(0) and ro < —^4(0) at h = (see table 
[II|). Tuning the value of f3, and starting from the same initial 
conditions, one can study both cases. 
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FIG. 3: Self-correlations for r > -A(0) in the T 
ferromagnetic quench. 



and h = 



In figures [3] and |4] we show autocorrelations for ro > 
— A(0) and ro < —^4(0) respectively. In both cases, the nu- 
merical estimations fit very well with the asymptotic predic- 
tions in Table HD In addition also the exponential factors char- 
acterizing the decay rates (r Q + A(0)) of these correlations 
are consistent with the analytical predictions. The difference 
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FIG. 4: Self-correlations for r < -A(0) in the T = and h = 
ferromagnetic quench. 

between the theoretical values for the parameters used and 
the numerical estimation obtained fitting the numerical data 
is smaller than 5 percent. 

At zero temperature, as we already mentioned, a system 
evolving from homogeneous ferromagnetic initial conditions 
is incapable to break the spatial homogeneity. In this con- 
text, the fact that autocorrelations decay to zero or non-zero 
constants according to the values of ro and A(Q) can be un- 
derstood in terms of the stability of the component <j){k = 0). 
For r > -A(0) one has from g3J that A(0) > 0. This 
in turn implies that the k = mode is unstable and the sys- 
tem evolves to a non-magnetized state. On the contrary, for 
ro < — A(Q) the stability condition .4(0) < drives the sys- 
tem to a non-zero magnetization. 

2. Zero Field Cooling 



The effects of the predicted critical field (lAlQft can be seen 
in figure [5] 




FIG. 5: Evolution of the magnetization patterns of a 100 x 100 por- 
tion of the 1960 x 1960 system for three different fields after a dis- 
ordered quench. First row: h = 0, second row: h = 4.0 < h c and 
third row: h = 4.7 > h c . Time grows from left to right. 



There, three dynamical sequences are shown after a dis- 
ordered quench for field values above and below the critical 
one. It can be seen that final configurations evolve, corre- 
spondingly, to modulated or homogeneous states as expected 
from our results (see Tabled and also predicted by thermody- 
namic equilibrium studies 1351 . 




200 400 600 800 1000 1200 1400 
At 



FIG. 6: Autocorrelations in the ZFC case for T — and h > h c 
and several waiting times. An asymptotic exponential relaxation is 
observed in agreement with analytic results in table U 

Autocorrelations for T = and h > h c decay very rapidly 
to a plateau which is dependent on t w . In figure[6]the quantity 
C eq (t w ) has been subtracted. It is seen that, as time grows, the 
exponential decays become evident, in agreement with analyt- 
ical results in table HI 




0.45 0.5 0.55 0.6 0.65 0.7 

(tt w ) 1/4 / ( t + t w ) 1/2 

FIG. 7: Collapse of the autocorrelations for the ZFC case, at T — 
and h = 0, using 12 samples. 

For h = 0, figure [7] shows a reasonable agreement with the 
predicted collapse of the two-time autocorrelations (see Table 
Hand [3 1]) in the asymptotic regime, attained for waiting times 
greater than 5 x 10 5 time steps. 

Finally, an important result concerns the growth law of do- 
mains at low temperatures, as discussed in the introduction. 
Our analytical results of eq. (TT3T) predict — t 1 / 2 , where 
£fc (t) is a positional correlation length. In figure[8]we show the 
time dependence of the azimuthally averaged structure factor, 
the positional correlation length and the orientational correla- 
tion length, defined in [13611. At long times all quantities seem 
to follow a logarithmic growth, at variance with the power law 
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with z = 2 found analytically. This slower regime may be due 
to pinning of topological defects, as observed also in a differ- 
ent model by Gomez et al. ll29ll . which is not captured by the 
Hartree approximation. More work is need to fully elucidate 
the origin of the logarithmic relaxation observed in figure [8] 




1000 10000 100000 1e+06 
t 

FIG. 8: Growing order in the disordered quench dynamics at T — 
and h — through the characteristic orientational and positional 
length scales (£nn(t) and £&(£)), and the peak of the azimuthal av- 
erage of the structure factor (S(km)), shown in a log-linear scale. 
Error bars are smaller than point size and the solid lines corresponds 
to logarithmic fits and are guidelines for the eye only. 



IV. CONCLUSIONS 

We studied the Langevin dynamics of a model with com- 
peting isotropic interactions at different scales at T = in 
the presence of an external field. Analytic results were found 
in the self-consistent field (Hartree) approximation for two- 
time autocorrelation functions and also for spatial correla- 
tions. Two different dynamical protocols (initial conditions) 
were considered: zero field cooled and field cooled. Analyt- 
ical results were confronted with numerical simulations of a 
prototypical system consisting of a ferromagnetic interaction 
frustrated with long range dipolar interactions in two dimen- 
sions. A detailed numerical scheme was developed and dis- 
cussed. 

Both analytically and numerically we found a critical field 
h c which divides the long time dynamic behavior of the model 
in two regions: for h < h c , the long time dynamics strongly 
depends on the particular initial condition, while for h > h c 
the initial conditions are irrelevant. The regime h < h c corre- 
sponds to the appearance of modulated structures at low tem- 
peratures (stripes or bubbles), while for h > h c the system is 
a paramagnet in an external field. 

We obtained explicit expressions for the autocorrelations 
and spatial correlation functions, both for the disordered and 
ferromagnetic quenches, considering all possible values of 
temperature and external applied field. In the FC case relax- 
ations are essentially exponential. In the ZFC case the system 
presents a slow coarsening of modulated domains. This coars- 
ening process reflects in an aging behavior of the two-time 
correlation functions. The analytical prediction for the aging 



scaling is in reasonable agreement with the numerical results. 
The results for the spatial correlations predict a growth law for 
the modulated domains with a power law £ ~ t 1 / 2 . This is not 
verified numerically, but a logarithmic growth £ ~ In t at long 
times is observed instead. This slow logarithmic behavior may 
be consequence of the very slow drift of topological defects, 
which are not taken into account by the self-consistent field 
approximation. 

Finally, it is important to remark that although the Hartree 
approximation does not consider the influence of topological 
defects in the dynamical behavior, our numerical results at 
T = compare quite well with the predictions of the theory. 
Surely topological defects are very important in the long time 
dynamics of the system and must be included in a complete 
theory of dynamics of modulated systems. Nevertheless, after 
a quench at zero temperature the dynamics of defects is ex- 
tremely slow, such that our approximate theory can still give 
valuable information on the relaxation in a time scale where 
defects can be considered as effectively "quenched". 
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Appendix A: Solution for K{t) at T = 

In this case equation © for K(t) takes the form: 

dK(t) 



dt 



= 2r K(t) + 2gV(t) 
+ 4ghK(t)^ 2 ^e- A{0)t x 

Jo 



2gh 2 



t i2 

dTKirf^e-^-^ 



(Al) 



Here a hypothesis is done as a very important step of the cal- 
culation that will allow us to linearize equation (1A1I) . namely: 

f K^^e-^-^dt' = (K 1/2 (t) 
Jo 

A(0)t > 1, (A2) 

and its consistency will be proved in the solution process, 
since the parameter ( is determined self-consistently using the 
hypothesis itself. By means of (IA2K equation (lAlb reads 

^ = 2r K(t) + 2gV(t) 

+ Agh foe-^ 1 (K(t) + 2gh 2 ( 2 K(t). (A3) 
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The term 4gh 0o e_A ^ ^ C^W * s unimportant in the long 
time limit, because it decays faster than others in the equation, 
given that A(0) > 0. Then, in the long time limit, we have to 
solve the following ordinary linear differential equation: 



dK(t) 
dt 



(2r + 2gh 2 ( 2 )K(t) + 2gV(t). (A4) 



The general solution of this equation is 

K(t) = K(0)e( 2ro+2 ^ 2 ^ 



+2g [ drV^e^ 2 ^^-^. (A5) 
Jo 

Different kind of solutions in the long time limit are possible 
according to the sign of 2ro + 2gh 2 ( 2 . 
If 2r + 2gh 2 ( 2 > 0, then: 



K(t) = ^(0) + 2^^°°dr V(r)e- (2ro+2 ^ 2c2)r 



x e v 



(A6) 



The integral converges because V(t) is a bounded, decreas- 
ing function of t. In this case we obtain a general solution 
independently of the particular function V(t). 
If2r + 2gh 2 ( 2 =0: 



K(t) = K(0) + 2g / drV(r). 



f' 

Jo 



(A7) 



This solution is less general than the previous one. It is valid 
while the term given by the higher order correction to the hy- 
pothesis (IA2I) is irrelevant in comparison with V(t). 

If 2ro + 2gh 2 ( 2 < 0, like before, there is not a general 
solution. The long time solution is given by the term decaying 
slower in expression (IA5K 

It is possible to obtain a critical field as the lowest applied 
field h c for which K(t) is an exponential. To do this we sup- 
pose that 2r + 2gh 2 ( 2 > 0. Then: 



K{t) = Ce^ 2 ^^\ 



(A8) 



and with condition (IA25 we obtain the following equation for 
the parameter £: 



c 



r + gh 2 C 2 + A(0) ' 



(A9) 



The limit case for which we have an exponential solution is 
when 2ro + 2gh 2 ( 2 = 0, applying this condition to the pre- 
vious equation we get ( c = 1/A(0). Then, the critical field is 
given by: 



h c 



ro 
9 



A(0). 



(A10) 



This field defines a zone in the parameter space (h > h c ) in 
which the long time behavior of the system is universal, in the 
sense that it is independent of initial conditions. 



1. h < h c 

We considered separately quenches from the disordered and 
ferromagnetic initial conditions. 



a. ZFC, disordered case 

With conditions ^ appropriately taken, the equation (IA1I) 
for the ZFC system reads: 



dK(t) 
dt 



2r K(t) + 2gAf(t) 



^l/2 e -A(0)(t-r) 



(All) 



+ 2gh 2 [ drK{r) 
Jo 

where, in the long time limit, f(t) is given by ll3Tll : 

f(t) = at~i, (A12) 



7 d - 1 „ d 

2tt2 / 7T 



with a = ^^y^-. 

We need to solve this equation considering the applied field 
to be lower than or equal to the critical one. In order to solve 
equation (lAlllh we use hypothesis (IA2I) obtaining the general 
solution: 

K(t) = K( ) e ( 2ro+2 ^ 2 ^ 

+ 2gaA f % e ^+2 9 h^){t-r) ^ (A13) 



Considering now 2ro + 2gh 2 ( 2 < (that is h < h c ) we get, 
in the long time limit: 



K(t) = 



-2gaA 



(2r + 2gh 2 ( 2 )t^ 

Applying condition (IA2b we obtain ( = , which allows 
us to rewrite K(t) as: 



K(t) = 



aAA 2 (0) 
{h 2 c - h 2 )ti 



(A14) 



When h = h c , higher order corrections to the hypothesis 
(lA2l ive: 



K(t) = 4gaA 



A(0) 



A(0) - 2r 



(A15) 



b. FC, ferromagnetic case 
Using the initial conditions ^ equation (lAlb becomes: 



dKjt) 
dt 



2 p -2A(0)t 



= 2roK(t) + 2g(t) 2 e 
+ 4^ e-^ f dr K^i^e-^^ 



2gh 2 



I drK(r) 1/2 e 
Jo 



l/2_-A(0)(t-r) 



(A16) 
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and, as before, using hypothesis (IA2I) it is possible to arrive to 
the following general solution in the long time limit: 

K(t) = K(0)e xt + 2g<g f dr e~ 2A ^ T e A( *- r) , (A17) 
Jo 

where A = 2r + 2gh 2 ( 2 . 

The analysis of the above solution when h < h c reveals 
different possibilities according to the different values of the 
parameters involved. In the long time limit we obtain: 

h = 



xt 



K(t) = Ce 



A is given by the solution of the equation 



(A18) 



A = 2r 



2gh 2 



K(t) 



(1 



r o +A(0) 

2g<\> 2 te- 2A ^ 



(L_) e 2r t forro > _^( ) 

for r = -A(0) 



(A(0) 



(A19) 



where hypothesis (IA2b has been used. Note that for all h < h c 
then A < 0, being A(ft c ) = 0. 



r o +A(0) 
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